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Abstract. The dynamical equation satisfied by the density matrix, when a quantum 
system is subjected to one or more constraints arising from conserved quantities, is 
derived. The resulting nonlinear motion of the density matrix has the property that 
the evolution is independent of the specific composition of the pure-state mixture 
generating the initial state of the system. 
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A nonlinear generalisation of quantum mechanics was proposed by Mielnik [U [2] , 
Kibble |3], and Weinberg [3] as an alternative to the linear evolution governed by the 
Schrodinger equation. The idea of this theory is to replace the linear Hamiltonian 
operator generating the evolution of the state vector with a state-dependent operator in 
such a way that the norm of the state vector remains constant. The properties of such 
"nonlinear observables" were subsequently studied in detail by Weinberg [5]. 

Following Weinberg's analysis, it was argued by Gisin |6l [7] and Polchinski 
[8] that the existence of dynamical nonlinearities in quantum mechanics might 
lead to undesirable physical features, such as the possibility of superluminal EPR 
communication. Indeed, it is often held now that the main issue associated with 
dynamics of the Mielnik-Kibble- Weinberg (MKW) type is that the evolution of a density 
matrix in general depends on the specific choice of pure-state mixture underlying the 
initial density matrix. Thus by observing the dynamics one might be able to determine 
the particular choice of mixture associated with the initial density matrix (Haag & 
Bannier [9]). 

Another potentially problematic aspect of nonlinear quantum mechanics was 
suggested by Peres pDj, who provided an example where the von Neumann entropy 
decreases in time. Tight experimental bounds on the deviation away from linear 
(unitary) evolution law have been found [11], while criticism of the MKW theory has 
been strengthened further by Mielnik [12] . As a consequence of these studies, one might 
conclude that a consensus has emerged to the effect that nonlinear quantum mechanics 
of the MKW type must be ruled out on physical grounds. 
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The purpose of the present paper is to report a generahsation of the MKW theory 
that circumvents some of these issues. The nonhnear dynamics that we propose emerges 
from the consideration of constrained quantum motions. We shall derive the dynamical 
equation satisfied by a density matrix when unitarity is compromised by the existence 
of one or more constraints. We demonstrate that although the resulting equations of 
motion are nonlinear, the associated evolution is autonomous and hence independent of 
the choice of initial mixture. We conclude that the nonlinearities arising from the type 
of constraints considered here can be regarded as representing a viable step towards an 
acceptable generalisation of the unitary evolution of quantum mechanics. 

Constrained motions appear not infrequently in the general study of dynamical 
systems [T3]. A systematic investigation of constrained motions in classical mechanics 
from a Hamiltonian point of view was carried out by Dirac [Ml [15]. In classical 
mechanics, the evolution is governed by a symplectic flow on phase space. The idea 
of Dirac, in essence, is to find the induced symplectic structure on the constraint surface 
in phase space, and to use this to characterise the dynamics. Dirac's approach has 
recently been applied in the quantum context to obtain the constrained dynamical 
equations satisfied by pure states [Ml H?! [T^ in various examples. In this paper we 
apply similar techniques to derive constrained equations of motion for mixed states, 
and to show that in the case of a pure-state density matrix the equations reduce to 
the results obtained previously. We consider in particular the special case where the 
constraints are given by the conservation of the expectation values of a family of mutually 
incompatible observables {^'^j (/c = 1, . . . , A^) such that [H, 7^ for all /c, where H 
is the Hamiltonian. Given an initial mixed-state density matrix po there are infinitely 
many different mixtures of pure states that give rise to po- In our scheme, nevertheless, 
the resulting nonlinear evolution for the density matrix is independent of the specific 
choice of mixture. We show also, for one version of the dynamics we consider, that the 
von Neumann entropy is a constant of the motion, and hence that the criticism of Peres 
does not apply to the constrained dynamics arising in that case. 

Let us now present our analysis. As remarked above, the constraints that arise most 
naturally in the case of constrained motions of density matrices involve the conservation 
of a family of observables {^^}k=i,...,N so that 

tr(p<l>'=) = (1) 

for all fc, where {c*^} are constants. For the moment we shall assume that N is even. 
The constraints are assumed to be non- redundant in the sense that [<l>'^, $'] 7^ for k 7^ / 
and [H^ 7^ for all k. We shall impose ([T]) by use of Lagrange multipliers, and write 

-^=i[p,^]-i5^A4p,<l1 (2) 

k=l 

for the proposed dynamics of the density matrix. Here {Afc}, k = 1,...,N, are the 
Lagrange multipliers conjugate to the constraints {c''}. In what follows we employ the 
usual convention that repeated indices are summed. 



N onlinearity and constrained quantum motion 3 

The object is to impose the constraints and derive exphcit formulae (given by 
nonlinear functionals of p) for the Lagrange multipliers. To derive the {A^} we consider 
the equation 



- tr(p$'^) = tr (^^ =0. (3) 

Substituting ([2]) in ([3]) gives 

tr([p, H]^'']) - A,tr([/9, l*^]*'^) = 0. (4) 

By the cyclic property of the trace operation, which says tr^ABC) = ti^CAB), this 
equation can be rewritten in the form 

tr(p[^, l''^]) - Ajtr(p[$^ l>*^]) = 0. (5) 

To proceed we define the antisymmetric matrix w^'^ by 

w^'' = tT{p[^^,^'']). (6) 

Provided that w^'' is nonsingular, which we assume holds at least initially, let us write 
Wij for its inverse so that Wijw^^ = 6^. Then the equation reads 

tr{p[H,^'']) = XjW^'' (7) 
and we can solve for the Lagrange multipliers to obtain 

Afc = ^/;,fctr(p[i^,<l>^]). (8) 
Substitution of this expression into ([2]) then gives 

^=i[p,H]-iw,Mp[H,^'mp,^% (9) 

This is the nonlinear equation of motion satisfied by the density matrix when it is subject 
to an even number of constraints of the form ([T]). 

Our objective next is to show that in the case of a pure-state density matrix the 
evolution equation iQ reduces in effect to the nonlinear Schrodinger equation obtained 
in reference [T7j. In particular, suppose that p is a time-dependent pure-state density 
matrix of the form 

for some state vector 1-0), not necessarily normalised. Then we see that 

and hence that 

^Pi i\ I ;\ mm I ;\ 



Therefore, for example, if p is assumed to satisfy the von Neumann equation 
dp 



(12) 



i[P,^], (13) 
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we can deduce that I?/)) satisfies the so-called projective Schrodinger equation 

\i^)-^^^\^) = -i{H-{H)m. (14) 
where {H) denotes the expectation of the Hamiltonian: 

The projective Schrodinger equation is essentially a slightly "weakened" form of the 
"full" Schrodinger equation 

|^) = -i^|^), (16) 

with the "unphysical" component of in the direction of removed. Clearly (|T6ll 
implies f|T^ . It is worth noting, however, that although f|T6|) is linear, the associated 
projective equation dH]), which embodies the physical content of the Schrodinger 
equation, is nonlinear. 

In the case of a constrained quantum system satisfying ([2]), essentially the same 
line of argument applies, and a short calculation shows that if p is a pure-state density 
matrix then 

with 



-i{H - {H)m - {^')m, (17) 



and where by analogy with ^ and ([H]) we have defined 



where tWjfc is the inverse of the matrix 



This result is indeed consistent with the equation derived in reference [TTj, and thus 
we are able to conclude that the dynamical equation ([2]) is a natural mixed-state 
generalisation of the constrained equation of motion for pure states. 

We remark, as we discuss in more detail below, that the motion generated by the 
constrained equation ([9]) depends only on the initial density matrix, and not on the 
choice of mixture leading to that matrix. This follows simply from the fact that ([9]) is 
autonomous in p. In addition, we can verify, by use of the cyclic property of the trace 
operation, that the von Neumann entropy, defined by 

S=-tr(plnp), (21) 
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is a constant of motion. The proof is as follows. First we note that if we let denote 
the components of p along the diagonal in a Hilbert space basis with respect to which 
p is diagonalised, we have 

dS* d 

-= -^tr(plnp) 

n 

= - ^PnlnPn 

n 

= -trf$lnpV (22) 



^dt 

Then by use of the dynamical equation for the density matrix we obtain 



^= - itr (^[i/-X;A,<l^p]lnpj 



= - itr \^\np,p] yH-J2>^k^' 

= (23) 

by virtue of the cyclic property of the trace, and the fact that p and In p commute. We 
can therefore regard Qj as representing a plausible candidate for an acceptable extension 
of the standard unitary dynamics defined by the von Neumann equation. 

We observe that the objections originally raised against nonlinear extensions of 
quantum mechanics were based implicitly on the essentially mistaken premise that the 
dynamics of a general mixed-state density matrix can and must be deduced from the 
dynamics of pure states. More precisely, if the initial density matrix /)(0) happens to 
take a decomposition of the form 

P(0) = $^p„fln(0), (24) 

n 

where {fln{.0)}n=i,... are normalised projection operators onto a set of pure states 
{\4'n)o}n=i,..., then it was assumed that the subsequent dynamics of the density matrix 
would have to be of the linear form 

PW = 5ZPnn„(t). (25) 

n 

The point is that in his original analysis Gisin ^ [7| had no way to deduce the dynamics 
of the density matrix except to regard it as following from the dynamics of pure states. 
From a modern perspective, however, we can take essentially the opposite view, and 
regard the density matrix as "fundamentally" representing the state of the system, from 
which properties of pure states can be deduced as special cases. Hence, in particular, 
there is no reason to suppose that the dynamics of p can or should be deduced, hnearly, 
from the dynamics of a set of hypothetical ensemble constituents. 
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So far we have considered the case for which the number of conserved observables 
is even. If is odd, then the antisymmetric matrix w^^ defined in ([6]) is singular and 
cannot be inverted. Thus, to obtain a system of constrained equations of motion that 
are apphcable to the case for which is odd, we need to modify the foregoing analysis. 
The idea here is to replace the commutator in ([21) by a symmetric product of the form 

^ = i[p, H] - A, ({p, ¥} - 2tr(p<l^)p) , (26) 

where {p, = p$'^ + denotes the anticommutator and where the {\k} comprise 
a set of Lagrange multipliers chosen to ensure the constraints ([T]), where need not be 
even. The plan is to circumvent the problem of the lack of invertibility of w^^ arising 
from the antisjTiimetric feature of the commutator by replacing it with a symmetric 
anticommutator {p, $'^}. The additional trace term on the right side of fl26l) is then to 
ensure conservation of the total probability, so that 

which follows at once from equation (!26|) . 

As before, we determine the Lagrange multipliers by considering the relation ([3]). 
Substituting fl26l) in ([3]), and using the cyclic property of the trace, we deduce that 

ii{p[H,¥]) = -iXj (^tr(p{<l>^ <!>*=}) -2tr(p<l^)tr(p<l'=)) . (28) 

To solve fl28l) for Xj we define the symmetric covariance matrix 

m^^ = tr(p{<|)^ <l^}) - 2tr(pl'^)tr(p<l'=). (29) 

Again, if m^^ is nonsingular, we can define its inverse ruij that satisfies mij-m^^ = 61- In 
this case, we obtain the following expression for the Lagrange multipliers: 

A,- = im.fetr (p[H, ^'']) (30) 



The invertibility of the matrix m-'^ is ensured by the positive-definiteness of the 
covariance matrix, which holds if p is itself nonsingular. Therefore, under this 
assumption the generalised equation of motion becomes 

^ = i[p, H] - im.fctr (p[^, ^']^ ({p, |.^}-2tr(p|.^)p) . (31) 

Example: As an illustration let us consider the case of a spin-| system for which 
H = dz and the single constraint observable is chosen to be $ = a^- The resulting 
motion for the density matrix is shown in Figure [H The space of density matrices in 
this case is a "Bloch ball" B of unit radius. The mixed states correspond to the interior 
points of B and the pure states form the surface of B. The chosen Hamiltonian would 
in the unconstrained case generate a rigid rotation around the z-axis. As we constrain 
the motion of the system by imposing the condition that the expectation of must be 
conserved, the resulting constraint surface corresponds to a cross-section of i3 at x = xq, 
where tii^pax) = Xq. From (13T!) we deduce the equation of motion for the system, and 
find, as is shown in Figure [1], that when the initial state is given by a pure state density 
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Figure 1. Constrained dynamics for mixed states in the case of a spin--^ system. The 
state space is a "Blocli ball" B of unit radius; pure states form the boundary (surface) 
of B and interior points represent mixed states. The Hamiltonian of the system is given 
hy H = &z- Hence under unitarity motion the evolution generates a rigid rotation of 
B around the 2;-axis. We constrain the motion by fixing the expectation of &x\ the 
resulting constraint surface is a slice of B perpendicular to the x-axis, as indicated 
above. Irrespective of the initial condition, the motion converges to the pure state 
for which ((t^) = asymptotically. The integral curves associated with the dynamical 
equation are shown in the figure. 



matrix, the state of the system remains pure as it evolves. In this case, the example 
reduces to the case considered in reference [19]. The mixed state evolution trajectories 
are also shown for a choice of initial states in the figure. The equator on the surface 
of B corresponds to a set of fixed states. Hence a state that initially lies at the point 
to the far left of the cross-section in Figure [T] remains fixed, and all other states evolve 
asymptotically towards the fixed point to the right in the figure. 

Returning now to (!3T]) . we note that in the special case for which p is a pure- 
state density-matrix of the form (fTOi) we are able to deduce the nonlinear projective 
Schrodinger equation satisfied by the state vector 1-0) that is applicable to both even 
and odd number of conserved observables. This is given by 

m - 1^1^) = - - ^^i^^' - (32) 

which is consistent with the result obtained in reference [19] . Therefore, (13T1) constitutes 
a natural generalisation of the result of [19] to the case of general density matrices. 

It is interesting to observe that, unlike the motion of (I9l), the motion of ( 13T1) . which is 
applicable to any number of constraints, does not necessarily preserve the von Neumann 
entropy. In particular, the entropy production is given by 

^ = -2AfeCov(<l>^lnp), (33) 
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where 

cov(X, Y) = tiipXY) - tr(pl)tr(pF). (34) 

The derivative of the entropy vanishes identically for pure states, which is why the pure- 
state limit (132|) is well defined. In general, however, we see that S is not necessarily 
constant. This is evident in the example shown in Figure [H On the other hand, just 
as in IQ, the evolution equation ( |3T1) is autonomous and independent of the specific 
composition of the mixture. Whether the fact that the entropy is variable raises an 

issue remains an open question. 
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